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ABSTRACT 

We develop new methods to study the properties of galaxy redshift surveys and radial 
peculiar velocity surveys, both individually and combined. We derive the Fisher infor- 
mation matrix for redshift surveys, including redshift distortions and stochastic bias. 
We find exact results for estimating the marginalised accuracy of a two-parameter 
measurement of the amplitude of galaxy clustering, A g , and the distortion parameter 
0. The Fisher matrix is also derived for a radial peculiar velocity survey and we discuss 
optimisation of these surveys for equal timescales. The Fisher Super-Matrix, combin- 
ing both surveys, is derived. We apply these results to investigate the 6 degree Field 
Galaxy Survey (6dFGS), currently underway on the UK Schmidt Telescope (UKST). 
The survey will consist of ~ 10 5 K-band selected galaxies with redshifts and a subset 
of ~ 15000 galaxies with radial peculiar velocities. We find for the redshift survey that 
we can measure the three parameters A g , V and (5 to about 3% accuracy, but will not 
be able to detect the baryon abundance, or the matter-galaxy correlation coefficient, 
r g . The peculiar velocity survey will jointly measure the velocity amplitude A v and 
r to around 25% accuracy. A conditional estimate of the amplitude A v alone, can be 
made to 5%. When the surveys are combined however, the major degeneracy between 
(3 and r g is lifted and we are able to measure A g , T, (3 and r g all to the 2% level, 
significantly improving on current estimates. Finally we consider scale dependence of 
r g and the biassing parameter b. We find that measurements for these averaged over 
logarithmic passbands can be constrained to the level of a few percent. 

Key words: Cosmology: theory - large-scale structure of Universe 



1 INTRODUCTION 

The quest to constrain cosmological parameters has been 
boosted in recent years by the emergence of large data sets 
from CMB experiments and galaxy surveys. In particular 
the first data release from WMAP combined with informa- 
tion from the 2 degree field (2dF) redshift survey and Lyman 
q forest data has allowed many cosmological parameters to 
be tightly constrained (Spergel et al. 2003). The Sloan Dig- 
ital Sky Survey (SDSS) when complete, will constitute a 
fourfold increase in redshifts over 2dF. The 6 degree Field 
(6dF) galaxy survey (Wakamatsu, 2003, Colless, 1999) offers 
a unique combination of a wide field redshift survey with a 
homogeneous subset of peculiar velocities. It is the aim of 
this paper, in the context of the present data-boom, to pre- 
dict the unique advantages of the 6dF galaxy survey over its 
contemporaries. 

Galaxy surveys have long been an invaluable source of 
information in cosmology, allowing cosmologists to infer the 
large-scale clustering of matter in the Universe, and cosmo- 



logical parameters. But despite their prominent role in our 
understanding of the clustering of matter, galaxy redshift 
surveys are fundamentally limited as a probe as we do not 
have a complete theory of galaxy formation from first princi- 
ples. This limitation is characterised by the galaxy bias pa- 
rameter, which is an unknown function relating the galaxy 
distribution to that of matter. In practice this function will 
be determined observationally, but its structure will help 
shape the theory of galaxy formation. 

Galaxy redshift surveys also suffer from the well-known 
distortion due to the use of Doppler redshifts to infer dis- 
tances. But redshifts are degenerate with peculiar velocities, 
the deviations from the Hubble flow generated by gravita- 
tional instability (Peebles, 1980). This projection of radial 
velocity information into redshift surveys provides aditional 
complication to the analysis of the clustering pattern, but 
also injects extra information on the dynamics of large-scale 
structure which couples directly to the mass distribution. 

In the local universe, the best way to probe the mat- 
ter distribution is from the local peculiar velocity field. At 



© 0000 RAS 



2 D. Burkey & A.N. Taylor 



higher redshift gravitational lensing becomes the method of 
choice for probing the distribution of matter, while on even 
higher redshifts the CMB provides the only viable and ac- 
curate method. 

By combining redshifts with distance indicators, such as 
the D„ — a relation, cosmologists have been able to estimate 
both true distance and radial peculiar velocities for local 
galaxies. Since estimating the true distance is a complicated, 
and time consuming process, these surveys have tended to be 
much smaller in size than the redshift surveys. In addition 
the distance estimators are subject to much larger uncer- 
tainties and biases than are the redshift estimates, and so 
have suffered from large uncertainties. 

One of the problems with this programme has been the 
difficulty in collecting a large homogeneous sample. Progress 
in this field has been limited by the need to patch together 
a number of surveys with almost inevitable systematic dif- 
ferences. The subject of the different biasses which plague 
peculiar velocity surveys is reviewed in Strauss and Willick 
(1995). There are two fundamental problems. D n -a and 
Tully Fisher relations derived from different data will be sys- 
tematically different, and a relation derived from one sample 
cannot be applied to another. Secondly, the observed quan- 
tities such as the velocity dispersions and apparent mag- 
nitudes needed for D n -a, will differ between surveys since 
each survey uses different observational methods and applies 
different corrections. The Mark III catalogue is one such 
compilation for which analysis has been done including the 
POTENT reconstruction where the 3-dimensional velocity 
is reconstructed from a potential - e.g. Kolatt et al. (2000). 
Another catalogue comprised of ~ 1600 field galaxies is the 
SFI catalogue. Both SFI and Mark III have been used to 
analyse the density power spectrum to constrain cosmolog- 
ical parameters as for example in Zehavi and Dekel (2000) 

In the age of mega-surveys this problem of compiling in- 
homogeneous surveys should be alleviated by large coherent 
galaxy samples with accurate distance indicators. The 6dF 
Galaxy Survey is the first such data set. The 6dFGS is cur- 
rently underway on the UK Schmidt Telescope (UKST). The 
survey will consist of ~ 10 5 K-band selected galaxies with 
redshifts and a subset of ~ 15000 of the brightest galax- 
ies with radial peculiar velocities. The galaxies are sampled 
from the K-band 2 Micron All Sky Survey (2MASS) Ex- 
tended Source Catalogue (Skrutskie, 2000) and so are dom- 
inated by early type galaxies. The advantage of the K-band 
is that it selects light from the old stellar population, and 
so presumably is a good indicator of the mass of the galaxy 
as a whole. 

In addition to redshifts the 6dF will also collect galaxy 
distances. The 6dFGS will be the first combined redshift 
and velocity survey, with the advantage that the selection 
criteria for both surveys are matched. The prospect of a 
homogeneous and integrated galaxy redshift and radial ve- 
locity survey is a big step forward in the analysis of peculiar 
velocities, which have suffered from both sampling and in- 
homogeneity effects. 

In this paper we develop an information-theory analysis 
of galaxy redshift surveys and radial peculiar velocity sur- 
veys, both individually and combined. Such an analysis is 
required both for survey design and to understand the sen- 
sitivity of the survey to cosmological parameters and so help 



determine which parameters the survey can be optimised to 
measure. These methods are general to the construction and 
analysis of galaxy redshift and radial velocity surveys, and 
can readily be applied to, e.g., the Sloan Digital Sky Survey 
(SDSS) dataset or the radial velocity fields probed by galaxy 
clusters from Sunyaev-Zel'dovich surveys of the CMB, such 
as in the case of the Planck survey (Lawrence and Lange 
1997). 

The paper is laid out as follows. In Section 2 we de- 
scribe our information theory methods, based on the Fisher 
information matrix. In Section 3 we derive the results we 
shall need for the analysis of galaxy redshift surveys while 
in Section 4 we apply these results to the problem of opti- 
mising their design. In Section 5 we estimate the sensitivity 
of a galaxy redshift survey to cosmological parameters, us- 
ing the 6dF as our fiducial model. The analytic results we 
require for an analysis of the radial peculiar velocity field 
are derived in Section 6, and used to optimise such surveys 
in Section 7. Parameter estimation from a peculiar velocity 
survey is studied in Section 8, while in Section 9 we analyse 
the combined surveys. We present our conclusions in Section 
10. There are also two appendices, with Appendix A deriv- 
ing a series of useful formulae for the analysis of redshift 
surveys, and Appendix B calculating the bivariate Fisher 
matrix. We begin with a review of the Fisher information 
matrix. 



2 FISHER INFORMATION ANALYSIS 
2.1 The single field case 

The Fisher information matrix measures the information 
content of a random field, 0(t"), about a set of parameters, 
0, that characterise that field. If this field is Gaussian dis- 
tributed, the Fisher matrix is (Vogeley and Szalay 1996; 
Tegmark, Taylor and Heavens 1997, Tegmark 1997; Taylor 
and Watts 2001), 

^ = \j 9 * 9j lnP^(fc) V eB (k) (1) 

where 

(0(fc)0*(fc)) = (2n) :i P t/>4 ,{k)8 D {k - k') (2) 

defines the power spectrum of a homogeneous field (f>, and 
we have assumed (<f)) = 0. The effective volume of the survey 
is 

where the spatial integral is taken over the volume sam- 
pling the 0-field. As the field is Gaussian all the informa- 
tion characterising the field is contained in the power spec- 
trum, P^ t j,(k) = P t f, r f,(k\d) . The gradients in equation (1), 
di = d/d6i, are taken in parameter space. 

If we randomly sample the field by a set of discrete 
points, the covariance of the sampled field is 

C^(k,r) = [P^(k)+N(r)] (4) 

where N(r) is a noise term that may vary spatially. 

The Fisher information matrix contains both the con- 
ditional error on a parameter on its diagonal and the 
marginalised error, 
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3 COSMOLOGICAL RANDOM FIELDS 



for the i th parameter. The Cramer-Rao inequality (Kendall 
and Stuart 1969), 



(6) 



where Ad = 6 — Go is the deviation of the parameters from 
their true value 0o, guarantees that the maginalised error 
is the smallest possible uncertainty on a parameter from a 
given measurement, and hence gives the minimum variance 
bound on marginalised parameter estimates. Furthermore 
the Cramer-Rao bound is an equality if maximum likelihood 
methods are applied. 

If the parameter is the power itself the uncertainty on a 
measurement of the power, band averaged over a logarithmic 
passband of width d\nk, is given by (Feldman, Kaiser and 
Peacock 1994) 



AP(jfc) 



27vP(k) 



^fk*d\nkV cS {k) 



(7) 



2.2 Multiple Fields 

Multiple fields can be incorporated into the Fisher Matrix 
formalism by generalising the definition of the data vector 
and its covariance matrix. Assuming we have two fields, (j>\ 
and (j>2 , we can form a covariance super- matrix 



C = 



Cll Cl2 

C21 C22 



(8) 



where each of the individual covariances C a p = (<j> a 4>p) 
where a, f) = 1,2 labels the fields. This super-matrix can 
be used to generate a Fisher super-matrix (eg Zaldarriaga, 
Spergel and Seljak 1997, Bouchet, Prunet and Sethi 1999), 



Tij= I ^fcdV^&CxfCovtCxCVF^CV. 



(9) 



The full expression of the Fisher super-matrix is given in the 
Appendix by equation (66) . 



2.3 Parameter correlation coefficient 

As well as the uncertainty on a measurement of a parameter 
set from a given survey, the Fisher matrix can also be used 
to estimate the correlations between measured parameters. 
The parameter correlation coefficient, 



la 



{A6,A6 3 ) 
AOiAO, 



Fr 1 



n 



(10) 



shows the degree of degeneracy between two parameters and 
allows one to identify the largest eigenvalue of the param- 
eter covariance matrix. This is of much use, as the degen- 
eracy between parameters can be the major limiting factor 
in parameter estimation. Knowing the principle degeneracy 
in parameter space, one can choose to increase the set of 
modes, preferentially adding information so as to decrease 
this degeneracy using the COMA data compression analysis 
discussed in Taylor et al. (2001). 



The design of galaxy redshift surveys has been investigated 
before by Heavens and Taylor (1997) who discussed the 
optimisation problem, and by Tegmark (1997) who calcu- 
lated the uncertainty on parameters in linear theory, ig- 
noring the effect on redshift-space distortions. Taylor and 
Watts (2001) generalised this to the nonlinear regime, and 
included both redshift distortions and nonlinear bias. Here 
we summarise these results for the linear analysis, assuming 
a general stochastic biasing model (Dekel and Lahav 1999) 
to relate the matter density field to the galaxy field. 

3.1 Galaxy density and radial velocity fields 

We shall assume for our analysis that the relevant field is 
the linear galaxy redshift-space density perturbation in the 
plane-parallel approximation. We write 



6 e g (k) = D(ka v n)[(b L + v?f)6 m (k) + e] 



(11) 



where / = din S/ din a w f2^ 6 (Peebles 1980) is the growth 
index of density perturbations, &l is a linear bias parameter, 
5 m (k) is the linear matter density field, fi = r.k is the cosine 
angle between the wave-vectors of the density field and the 
observers line of sight and e is a stochastic bias term. We 
assume (e) = and (e 2 ) = a 2 . We have also included a 
nonlinear redshift distortion term 



D{ka v n) = (l + k 2 a 2 v n 2 /2y 1/2 



(12) 



where a v is the pairwise radial velocity dispersion. 

We define a second bias parameter as the ratio of power 
in the galaxy distribution P gg (k), to that in the matter 
fields,P mm (k) : 



b = 



P 99 {k) 



bl + 



(/,•) ! /> C (kY ^ 

We may also introduce a galaxy correlation coefficient, r g , 
defined by 



Pgmik) 



\J Pgg (^)Pmm (&) 



b 



(14) 



(Dekel and Lahav 1999) which will allow us to change be- 
tween the two bias parameters. Although &_t is more physi- 
cal, it is b which is more commonly measured. 

We choose our second field to be the radial gradient of 
the radial peculiar velocity field, defined in Fourier space as 



u'(k) = -^-f.v(k) 
or 



- M H/(n m )5 m (fe) 



(15) 



where H is the Hubble parameter. The modes of the ra- 
dial velocity are uncorrelated with the modes of the density 
field. However the Fourier space radial gradient of the radial 
velocity is correlated and will yield a cross power spectrum. 
The auto- and cross-power spectra of these two fields 

are 



Pggik) = 

Pu'u 1 (fe) 
Pgu'(k) 

where 



D 2 {ka vjJ ){\ + 2 f i 2 r g (3 + //3 2 )6 2 P mm (fc)(16) 

= fl 4 H 2 f3 2 b 2 P mm 

= -ii 2 H(3D{ka v ii){r g +^(3)b 2 P mm {k), (18) 



(*). 



(17) 



(19) 
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is the linear redshift distortion parameter. 

The three spectra to be measured are the redshift-space 
galaxy power spectrum, P gg (k), the radial velocity gradient 
power spectra, P u i u i(k), and the cross-spectra of these fields, 

The noise terms associated with these fields are 

N 99 (r) = ^r y (20) 

N u , u ,{r) = ^ all{r) 



n v (r) 



N gu ,(r) = 0, 



(21) 
(22) 



where n g (r) and n v (r) are the number densities of the galaxy 
and velocity surveys, respectively. The factors of /i 2 and k 2 
are due to our use of the radial velocity gradient. We quan- 
tify the number densities in Sections 4 and 7. The noise term 
for radial velocities arises from the intrinsic uncertainty on 
galaxy positions due to the dispersion in the D n -a relation 
for ellipticals. This can be approximated by 

a DI (r) = a Hre r/r °", (23) 

where ao is the fractional error on the distance indicator. 
We allow the fractional distance error to blow up beyond 
some fiducial distance, r orr . We choose a conservative r crr of 
135/i _1 Mpc (which is 90 % of the 6dF median depth), but 
allowing this to decrease to as low as 50% of the median 
depth had little effect on our results. 



3.2 The matter power spectrum and cosmological 
parameters 

Throughout we define the matter power spectrum by 

Pmm(k) = Q 2 k n T 2 {k;V,uj b ), (24) 

where Q is the amplitude of matter perturbations, n is the 
primordial spectral index of perturbations and T(k) is the 
matter transfer function. The dependent parameters of the 
transfer function are the so-called "shape parameter" , 

T = fl m h (25) 

which stretches the scale of the transfer function, and 



(26) 



where f2& is the density parameter of baryonic matter. Here 
we use the form of the transfer function given by Eisenstein 
and Hu (1998). 

The other parameters which will be of interest to us are 
the amplitude of galaxy clustering and the amplitude of the 
velocity field. As we shall assume linear theory throughout, 
these are respectively 



A g = bQ 



(27) 



where 6 is the linear galaxy bias parameter, and 

A v = f(n m )Q- (28) 

Also of interest is the mass-galaxy correlation coefficient r g , 
and as we shall assume the galaxy survey is in redshift- 
space, the final parameter is the redshift-space distortion 
parameter, /3. 

For the purposes of this paper we shall assume the 
primordial spectrum is scale- invariant so that n = 1. 



Hence the parameters we shall investigate are the six non- 
degenerate parameter set (A g , A v , T, uib, f3, r g ). The model 
we shall use has fiducial values of (Q,h,£l m ,t+>b,b,r g ) = 
(5 x 10~ 5 , 0.65, 0.3, 0.025, 1, 1) giving the measurable param- 
eters the values (A g , A v , T, u>b, (3, r g ) = (5 x 10~ 5 , 0.486 x 5 x 
10~ 5 , 0.195, 0.025, 0.486, 1). 



4 THE GALAXY REDSHIFT SURVEY 
4.1 Survey parameters 

A galaxy redshift survey can be parameterised by the frac- 
tional sky coverage, /sky, and the mean radial density of 
galaxies in the sample, which we parameterise by 



n g (r) = N g 



rlT{3/a] 



-(r/r,) a 



(29) 



where r* is a convenient survey depth scale, N g is the total 
number of galaxies and typically a — 1.5. 

Rather than use r* we use the median depth of the 
survey (r m m 1.5r») and the number of galaxies available in 
the survey roughly scales as 



AT, « 4. x 10~ 2 ar2; 125 , 



(30) 



from fitting the near infrared differential number counts 
compiled in Kochaneck et al. (2000). 

The sampling fraction, a, is simply the fraction of 
sources expected at any redshift which are actually observed. 
We shall assume that space is Euclidean. The three survey 
parameters are then /sky, r m and a. 

4.2 Optimisation of a Galaxy Redshift Survey 

The problem of optimisation of a galaxy redshift survey has 
been previously considered by Heavens and Taylor (1997), 
who showed that the optimisation for a fixed timescale re- 
duces to a one-parameter optimisation problem. We first 
need to choose a parameter for which to optimise. Here we 
choose the amplitude of galaxy perturbations, A g , where the 
conditional error on a measurement is given by 



AA a 



1 f_£i 

2 J (27T 



-1/2 



(31) 



We now assume the timescale for the redshift survey is 
a constant and scales as 



t = t af sky (r m /h ^pc) 7 , 



(32) 



where to is some fixed timescale. We take as our baseline 
a 2 year survey where t/t ~ 7 x 10 17 . The strongest con- 
straint then comes from the survey depth. The uncertainty 
on parameters scales as A8 oc f sk J 2 , while the time scales 
as t oc /sky, indicating, as is well known, that as large a solid 
angle as possible should be chosen. With this time constraint 
the survey median depth is 



r m = 150(a/sky)" 1/7 ^ 1 Mpc. 



(33) 



This states that the depth of the survey is increased at the 
expense of sparse sampling or reduced sky coverage. We can 
now minimise the uncertainty on A g . 

Figure 1 shows the results of minimising the fractional 
error on the galaxy clustering amplitude as a function of 
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Figure 1. Optimisation of the 6dF redshift survey. Here we plot 
the fractional error on a conditional measurement of the ampli- 
tude of galaxy clustering, A g , as a function of survey median 
depth, r m . The lines are for constant survey timescales with sky 
coverage of unity (bottom), three quarters (second from bottom), 
a hemisphere (third from bottom) and a quarter (top). 

survey depth, for a fixed timescale. Each of the lines is for 
a different sky coverage. The approximate relation we use 
between limiting K band magnitude and median depth is 

KumU * 5los (iFm^) +L865 - (34) 

The sampling fraction for each value of r m may be calculated 
from equation (33). The dashed line in Figure 1 shows the 
physical limit a = 1. The optimal survey is a shallow all-sky 
survey with depth r m = 140/i _1 Mpc. For a hemisphere this 
increases slightly to r m = 150/i~ 1 Mpc, while the optimal 
value of a ~ 0.7 is fairly insensitive to either f 3 k y or r m . The 
increase in the uncertainty in the fractional error in A g with 
increased depth is due to the decrease in the total number 
of galaxies in the survey. At the extreme of very large r the 
fractional uncertainty in A g tends towards A g ~ r 6,5 . 

Figure 2 shows the linear galaxy power spectrum using 
the optimal redshift survey parameters for the 6dF, along 
with the expected conditional error bars in band averaged 
pass-bands of width A In A; = 0.5. We have also plotted the 
effective volume, V® g (k), which indicates where the informa- 
tion is maximised in the survey. Finally we have also plotted 
the noise per mode, k 3 /2ir 2 n, where n is the average density 
of the survey. For the rest of the analysis we use Ku mit = 12, 
fsky = 0.5, r m = 150/i _1 Mpc and a — 0.7. 



5 PARAMETER FORECASTS FOR THE 6dF 
REDSHIFT SURVEY 

In linear theory a galaxy redshift survey based on a CDM 
model contains three free parameters: the amplitude of the 
galaxy power spectrum, A g — bQ, the scale of the break, 



q [ i i i i i i i I I I I I 

no -3 o.oi o.i l 10 ioo 

k/h Mpc -1 

Figure 2. The linear galaxy power spectrum with fiducial LCDM 
parameters (solid light line). Also plotted are the expected band 
averaged error on measured power pass-bands (grey area), and 
the effective survey volume V^ s (k) (dark solid line) and the shot- 
noise per mode (dashed line), for the 6dF redshift survey. We 
parameterise the 6dF survey by its sky coverage, depth and total 
number of galaxies. 

parameterised by V — Q. m h and the redshift distortion pa- 
rameter, /3 = S7°' 6 /&. Other parameters that could be con- 
sidered are the small-scale pairwise velocity dispersion, a v , 
or the baryonic density parameter, fib, which can introduce 
oscillations in the the matter power spectrum. Other possi- 
ble parameters include the spectral index n (Tegmark 1997) 
or the correlation between luminous and dark matter r g . We 
consider the estimation of five parameters A g , V, j3, fib and 
r s- 

5.1 The information content of the redshift-space 
galaxy power spectrum 

The information content of a galaxy redshift power spec- 
trum can be examined per Fourier mode by plotting the 
derivatives of the log-power with respect to the parame- 
ters, din P gg (k) / d6. In geometric terms this derivative can 
be thought of as a vector in Hilbert space, and the Fisher 
matrix as the dot product of such vectors. Hence the more 
similar two curves are, the more "parallel" they are and so 
the greater the degeneracy between the parameters. In the 
extreme that two curves are identical the Fisher matrix be- 
comes singular and the uncertainty becomes infinite. 

Figure 3 shows that the derivatives for the five parame- 
ters we can extract from linear theory, 6 — (A g , F, /3, u>h,r g ). 
The Ag-derivative is a constant across all modes, as ex- 
pected. Similarly both r g and (5 are almost constants, as 
they appear in the normalisation of the galaxy power spec- 
trum. Hence we can expect some degeneracy between these 
three parameters. 

The final two parameter derivatives, F and u>b are also 
similar. This is because F parameterises the break in the 
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Figure 3. Derivatives of the matter power spectrum 



CDM power spectrum, between retarded modes inside the 
horizon at matter-radiation equality and growing modes out- 
side, while Ui parameterises the effect of baryonic damping 
of perturbations inside the horizon when the baryonic fluid is 
oscillating. Hence both parameters are related to two phys- 
ically different effects that produce similar damping of the 
matter power spectra. An effective shape parameter can be 
defined by (Peacock and Dodds 1994) 

r cff = Ye- 2a " h (35) 

which crudely models both effects but explicitly shows the 
degeneracy. On scales larger than the break, the effects of 
the baryons dies away quicker, and so on the largest scales 
the power spectra contains more information on F than on 

U>b- 



5.2 The amplitude of galaxy clustering and 
redshift distortion parameter 

In the limit of negligible shot noise, it is possible to calculate 
analytically the two-parameter Fisher matrix for A g and 0. 
The details of this are shown in Appendix A. From equa- 
tion (54) the correlation coefficient between the two param- 
eters is "iAgp = —0.78, which is very close the value found 
by Taylor et al (2001) for the Point Source Redshift Survey 
(PSCz). The strong anti-correlation is due to the similar 
effects of the parameters on the power spectrum, as they 
both appear in the amplitude of the redshift-space power. 
From equation (52) we can read off the fractional errors: 
AA g /Ag = 0.01, and A/3//3 = 0.057 for current typical sur- 
vey sizes at the limit of linear theory k w 0.2/iMpc . 

For a single parameter estimate of j3, again in the ab- 
sence of shot noise we find; 



Ap 

P 



2tt 2 



3tan- 1 ( A /^) 



-1/2 



(36) 



1+/3 VP 
where V is the survey volume, and k is the maximum 




q 



d 



d 

I I- 



I 1 1 1 1 1 



W g 



H 1 — I I I I | 



?(A g r) 

7(A g 0) 

7(r/s) 



o — 



b.o 



1 



0.02 



0.05 



0.1 



0.2 



k/h Mpc 1 

Figure 4. Upper plot: The expected fractional uncertainty on 
A g , r and (3 from the 6dF redshift survey. Lower plot: The cor- 
relations between A g , V and f3. We model the survey as described 
in the text. 



wavenumber used to estimate 0. For (3 = 0.5, and k 
0.2 h~ Mpc this gives a fractional uncertainty of 

-1/2 



A? 



V 



4.5 x 10 4 [/i- 1 Mpc] 



(37) 



Hence, even for a well-sampled galaxy redshift survey, the 
effective volume has to be large to beat down the sam- 
pling noise on a determination of (3. To achieve a 1% error, 
assuming all other parameters are known, requires a sur- 
vey with volume V ~ 4.5 x 10 8 [/i~ 1 Mpc] 3 , over four times 
larger than the 2dF, 6dF or SDSS redshift surveys. Hence 
to achieve high accuracy on cosmological parameters in the 
linear regime requires an order of magnitude increase in sur- 
vey size. However even for these surveys, with volumes of 
order V « 10 8 [/t -1 Mpc] 3 , we can expect an accuracy of 
A/3/f3 ~ 0.02. An alternative to going large is to extend the 
analysis into the nonlinear regime, as outlined by Taylor & 
Watts (2001). Finally, the uncertainty can be reduced by 
combining data sets, as investigated here. 
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5.3 Three parameter set: A g , V and f3 

Introducing extra correlated parameters into the analysis 
causes uncertainties to increase, and adds extra complex- 
ity to the analysis. Hence we first consider a reduced 3- 
parameter set consisting of A g , F and /3, and then see how 
adding new parameters changes our results. 

Figure 4 shows the expected fractional marginalised un- 
certainties from the 6dF redshift survey plotted as a function 
of maximum wavenumber analyzed. The scale of interest for 
linear analysis is around k ~ 0.2 /iMpc -1 . At this wavelength 
the best estimated parameter is V which is constrained to 
within about 2%. The uncertainty on A g and f3 is about 3%. 
For f3 this is close to our conditional estimate in Section 5.2. 

Figure 4 also shows the correlation coefficient for these 
parameters. The strongest correlation is between the am- 
plitude of galaxy clustering, A g , and the shape parameter, 
r, with 7(j4 s r) ~ —0.95. This arises because a change in 
amplitude can be mimiced by a shift in scale. The least cor- 
related parameters are F and /3. 

Intermediate is the correlation between A g and ft. This 
is of particular interest as their combination gives (Taylor 
et al 2001) 



A g p = Qf(n r , 



(38) 



and so provides an independent estimate of the amplitude 
of matter perturbations. The correlation of these parame- 
ters can be calculated in the limit of negligible shot-noise, 
and ignoring other parameters is given by equation (54) in 
Appendix A. For /3 = 0.5 we find 7a,/3 = —0.78, very close 
to the value we find for a two parameter analysis. 

As these parameters are anti-correlated, the combina- 
tion A g (3 has a smaller formal error than either A g or (3. 
This arises because the combination A g /3 marginalises over 
the A g -/3 parameter space along the longest line of the de- 
generacy, so that the marginalised combination has a smaller 
error than its components. This effect is shown later, in Sec- 
tion 9 figure 11. 



5.4 Four parameter set: A g , F, /3 and u>t 



Introducing a fourth parameter, uib, has little effect on the 
uncertainties on A g and [3 on scales k > 0.05 /iMpc" 1 , as 
shown in Figure 5. Figure 5 shows that this is due to the 
weak correlation between u>b and these parameters. How- 
ever the strong correlation between ujb and T, ^/(Tujb) ~ 0.9, 
which we expect from our discussion in Section 5.1, degrades 
the uncertainty on F. Disappointingly the uncertainty on a 
measurement of cot is large, Aujb/uJb ~ 0.45, which reflects 
its weak effect on the matter power spectrum. As with all 
the parameters, the main way to improve on sensitivity is 
with a larger redshift survey. 

Since ujb and Y are correlated would we do any better 
by estimating the baryon fraction, 



(39) 



rather than simply cj(,? Unfortunately not. Although these 
parameters are positively correlated, and so the marginalised 
combination will project out along the longest axis of the 
u)b-T error ellipse, making the fractional error on fb smaller, 
other correlations, particularly between u>b and A g , prevent 
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Figure 5. Upper plot: Predicted errors for a four parameter 
analysis involving u> b . Lower plot: Correlations in a four param- 
eter space with u>b 



any improvement. We have re-analysed the survey with ft 
replacing Ub throughout and find negligible change in our 
results. We conclude that the main limiting factor here is 
the survey volume. 

In order to reduce the error on LUb to the 10% level we 
predict one would require a survey some 10 times larger. 



5.5 Four parameter set: A g , F, f3 and r g 

Instead of uj b , we could have chosen r g as our fourth param- 
eter. Figure 6 shows that this time A g and T are unaffected 
by this new parameter, but the uncertainty on f3 is dramat- 
ically increased, with both the fractional errors on (3 and 
r g around 35%. Again this can be traced to the expected 
degeneracy between r g and (5 discussed in Section 5.1, and 
shown in Figure 6, where 7/3,- g ~ — 1. Interestingly the cor- 
relation between A g and /3 has now become positive, so that 
the error on their combination will be increased. 
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Figure 6. Lower plot: Uncertainties for a four parameter anal- 
ysis invloving r g . Lower plot: The correlations between A g , T, 
/3 and r g 



6 GALAXY VELOCITY SURVEY 

The application of Fisher information methods to galaxy 
velocity surveys has so far been limited. Here we derived 
for the first time the Fisher matrix from parameters from 
such a radial velocity survey, and use it to find the opti- 
mal survey parameters, estimate the conditional errors on a 
measurement of the velocity power spectrum and estimate 
the marginalised errors on cosmological parameters. In Sec- 
tion 9 we combine the radial velocity and galaxy redshift 
surveys. 

Figure 7 shows the 3-D velocity power spectrum for our 
fiducial LCDM cosmology, and the effective survey volume, 
k 3 V t ^{k)/2n 2 , for a model 6dF velocity survey. 
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„-^ 3 Y eff (k)/2TT £ 
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no ° 0.01 0.1 1 10 

k/h Mpc -1 

Figure 7. The peculiar velocity power spectrum, A v (k) = 
k 3 P v (k)/2n 2 with predicted noise and effective volume for the 
6dF peculiar velocity survey. The velocity survey is parameterised 
by /sfej/i cut-off depth, a v and the accuracy of the distance indi- 
cator. 

7 VELOCITY SURVEY OPTIMISATION 

7.1 Velocity survey parameters 

The velocity survey can be parameterised by the fractional 
sky coverage, / s k y and the mean radial density of galaxies 
in the sample, which is parameterised as for the galaxies 
(equation 29). The early type galaxies for the velocity survey 
will be a subset of this group and so for analysis of the 
velocity survey the number count needs to be multiplied by 
a factor of 0.3 and the analysis limited to r m which is the 
cut-off in the 6dF velocity survey. We shall again assume a 
Euclidean universe for simplicity, and as we are working at 
low redshift. The survey can then be parameterised by / s k y , 
r m , a v and ctq, the expected fractional error in the distance 
indicator. 

The angular part of the effective volume for the velocity 
survey can be integrated analytically and yields 



V3t(k) = 2nf aky / drr 2 Q 



n u (r)P uu (k) 



(40) 



where the function O(z) is defined in Appendix A, equation 
(50). Figure 7 shows the dimensionless effective survey area. 
The peak of the effective volume again indicates where the 
information content is maximised. 



7.2 Band-averaged power 

As well as showing the effective volume Figure 7 shows the 
error bars expected for a measurement of band-averaged ve- 
locity power in a series of pass-bands, using equation (7). 
The power spectrum plotted is the three dimensional ve- 
locity power spectrum but the error bars reflect the fact 
that this will be derived using only radial velocity data. The 
pass bands are of width A In k = 0.5 and so this plot may 
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be directly compared with the corresponding band averaged 
power for the redshift survey in Figure 2. The velocity power 
spectrum has much larger error bars which is a reflection of 
the smaller size of the velocity data set. 

7.3 Optimisation of a velocity survey 

As with the redshift survey we wish to know whether the 
velocity survey can be optimised in the sense of maximising 
the information content for the minimum observing time. To 
optimise the survey, we choose the amplitude of the power 
spectrum as a fiducial optimisation parameter, where 

1 -1/2 

(41) 



AA V 

A v 



din k 



4tj-2 



and we consider the case of equal-time surveys. 

The timescale for a velocity survey scales as the total 
number of sources in the survey, N ~ r 3 , the reciprocal of 
the effective flux of the sources, S ~ r -2 , the distance error 
per source, odi, the fraction of the sky surveyed, /sky, and 
the fraction of the sources sampled, a. The timesale for a 
velocity survey then scales as 



t » Ovft 



sky 



f_£!»_y 



h~ 1 Mpc 



to 



(42) 



where to is a constant timescale. For a survey of a few years 
we find the ratio t/t ~2x 10 13 . 

We can now minimise the fractional uncertainty on the 
amplitude of the velocity power spectrum, A v , in survey- 
space with the equal-time constraint and the restrictions 
/sky, ctv < 1. In our model, the overall uncertainty on param- 
eters scales as A8 oc / sky , so all our results and timescales 
can be scaled by the fraction of sky surveyed. In the following 
analysis we fix the value of / sky and optimise in terms of the 
remaining survey parameters. In addition the uncertainty 
in source distance per source, oo, and the source sampling 
rate, ot v , always appear in the ratio ag/a v , further reducing 
the number of free survey parameters we need to consider. 
Finally, we can relate the remaining parameter, r m to the 
other parameters by 

CO 



150 



sky 



Vkms-V 



h 1 Mpc. 



(43) 



Hence we only need to optimise the survey with respect to 
the degenerate parameter 



1 I to \ 
JaZ Vkms -1 / 



(44) 



In practice there is a lower limit to the accuracy of the dis- 
tance indicator. The current best accuracies from the D n -a 
relation provide distance estimates with a minimum scatter 
of (T m i n ~ O.lkms -1 e.g. (Jorgensen et al. 1993; Gregg 1995), 
and so we shall optimise the survey design in terms of the 
new parameter 



Tj - 



1 / (TQ + a m in \ 

faZ V kms -1 / 



(45) 



which better reflects the survey's limitations. 

Figure 8 shows the fractional uncertainty on a condi- 
tional measurement of A v for equal-time surveys as a func- 
tion of r) and where the different curves correspond to differ- 
ent sky fractions, / sky . The lowest, solid line is for an all-sky 




7 7=( CT o+ CT min)/« 1/2 (kms- 1 )" 1 

Figure 8. The normalised fractional uncertainty on the ampli- 
tude of the velocity power spectrum for the degenerate parameter 



,) 1 



for equal-time surveys. The four lines corre- 



kms 

spond to different sky fractions (see Figure 1). 



survey, with higher lines reduced by one quarter. The ex- 
pected 6dF velocity survey lies somewhere between the top 
two lines depending on dust extinction from the plane of the 
galaxy. 

The other parameter that this analysis will depend on 
is the upper wavenumber used in the analysis, fc max . Lin- 
ear theory is reasonably robust up to k — 0.2 h~ Mpc. In 
practice it was found that the optimal survey parameters 
changed little with changes in k max , so that our results with 
kmax = 0.2/tMpc -1 are quite general. 

Our results show that the optimal strategy is to be as 
accurate as possible in the distance determinations at the 
expense of going deep. In fact the absolute optimal strat- 
egy is found to be impossible with current distance indi- 
cator limits - if the distances could be better constrained 
then we could do even better. The reason for this is that by 
k « 0.2/iMpc" 1 the power spectrum is already becoming 
dominated by shot noise. If there were more sources (and so 
less shot noise) this would have the same effect as increasing 
a in the degenerate parameter rj and so the optimal value 
for Co would also increase. The bottom line is that the sur- 
vey should be as well sampled and as accurate as possible 
before going deep. 

For the rest of this analysis we assume a — O.lfcms -1 , 
a v — 0.2, f sky — 0.5 and r\ — 0.2. 



8 COSMOLOGICAL PARAMETER 

FORECASTS FROM VELOCITY SURVEYS 

In the linear regime the properties of the velocity field can 
be specified by just two parameters - the amplitude of the 
velocity power spectrum, A v , and the shape parameter, V. 
Figure 9 shows the fractional error expected on A v and T for 
an optimal survey, as a function of maximum wavenumber 
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Figure 9. The expected fractional error (top) and correlation 
coefficient (bottom) for the amplitude and shape parameter of 
the 6dF velocity power spectrum. The survey parameters are dis- 
cussed in the text. 



k/h Mpc 1 

Figure 10. Supcrmatrix fractional errors and correlations be- 
tween A g , r and j3. These predictions are for a combination of 
the 6dF redshift and peculiar velocity surveys. 



analyzed. If we truncate the analysis at k = 0.2 ft. - Mpc, 
the fractional uncertainty in both A v and F is expected to 
be ~ 25%, ignoring any systematic uncertainties. 

Figure 9 also shows the parameter covariance coefficient 



7A„r(fc) 



(AA V AF) 



-0.8 



(46) 



AA V AV 

as a function of wavenumber. It is because of this anticor- 
relation that a joint constaint of A v and V is difficult. If 
we assume a value for V obtained from the redshift survey, 
Figure 8 shows that A v may be constained on its own to 
« 5%. 



9 COMBINING DATA SETS 

In this section we combine the information content of the 
radial velocity field and the galaxy redshift survey. We cal- 
culate the Fisher matrix using equation (9) which is written 
explicitly in Appendix B. As explained in Section 3 we now 



work with the radial velocity gradient as this will yield a 
cross-power spectrum with the density field. 



9.1 Three parameter set 

To see how combining redshift and velocity data affects pa- 
rameter estimation we begin with our three parameter set; 
A g , r and /3. Figure 10 shows the uncertainties and correla- 
tions of this set. Comparing with just redshift information, 
in Figure 4 we see that the constraints on A g and F remain 
largely unchanged and there is a slight improvement in the 
constraint on f3. So for this three parameter set there is rel- 
atively little improvement. This is because the inclusion of 
the velocity data does nothing to break the main degener- 
acy between A g and T. The other parameters were relatively 
uncorrelated. The lower panel in Figure 10 shows that the 
correlations remain unchanged from the redshifts-only anal- 
ysis. 

Figure 11 shows the A In (Likelihood) = — | likelihood 
contours about the fiducial maximum likelihood point in a 
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Figure 11. lcr contours for likelihood marginalised over one free 
parameter. The outer contour is for the redshift survey only. 
Hatched regions indicate the constraint from combining with the 
peculiar velocity survey. 

two-parameter space. These contours represent the likeli- 
hood marginalised in each case over the one remaining pa- 
rameter. The outer contour is that from the redshift sur- 
vey alone whereas the inner hatched region is the constraint 
when velocity data is included. The maximum wavenumber 
analysed here is k — 0.2 ft _1 Mpc. 

9.2 Four parameter set 

As with the galaxy survey, we can add in a new parameter 
and study its effect on a likelihood analysis. Adding in u)b 
has little effect on A g and /?, but degrades F, as was the 
case for the redshift survey alone. The baryonic density is 
not detected from the combined survey anyway, nor is the 
baryon fraction. We shall not consider this set further. 

Figure 12 shows the marginalised uncertainties of a four 
parameter set including the galaxy-mass correlation param- 
eter, r g . Comparing with Figure 10 we see that A g and Y 
are unchanged by the inclusion of the new parameter, since 
these are practically independent of r g . The redshift distor- 
tion parameter (3, has degraded slightly, due to its strong 
correlation with r g , but is still measurable at the 3% level. 
The big change from the redshifts-only analysis is in the 
galaxy-mass correlation parameter r g , which is now measur- 
able from the combined surveys, with a formal uncertainty 
of better than 2%. Comparing Figure 12 with its redshift- 
survey-only equivalent Figure 6 it is clear that there has 
been a great improvement in the joint constraint of and 
r g . The reason for this can be seen by comparing the cor- 
relations between the two parameters in the lower panels 
of both figures. In the redshifts-only analysis (5 and r g are 
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Figure 12. Supermatrix fractional errors and correlations be- 
tween A g , r, (3 and r g . 

strongly anticorrelated whereas when the redshifts are com- 
bined with peculiar velocities this degeneracy is broken and 
the parameters become almost independent. 

Again we have plotted the expected likelihood contours 
for this analysis in order to help visualise the improve- 
ment in the constraint. In Figure 13 the contours represent 
A In (Likelihood) = — | after marginalization over the two 
remaining parameters. Again the hatched region shows the 
improvement from combining the data sets. Particularly in- 
teresting is the panel showing contours for (3 and r g which 
illustrates the breaking of the degeneracy between the two, 
and the subsequent improvement in their joint constraint. 

9.3 Scale dependence of b and r g 

Given the tight constraint on r g and j3 it is interesting to 
consider a possible scale dependence in either r g or the bi- 
assing parameter b. Although we have approximated b to be 
constant over the scale of interest, it may have some scale 
dependence. Another way to consider scale dependence in 
the biassing is to fix 6 as constant and allow the parameter 
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Figure 13. lcr contours for likelihood marginalised over two free 
parameters. 



r g to vary with k. The upper panel of Figure (14) shows the 
expected error bars on a measure of r g band averaged over 
a logarithmic passband of width Alnfc. The width Alnfc is 
marked on the plot as a series of dots, b is fixed at b — 1 for 
this analysis. The results are extremely encouraging; r g (k) 
can be measured to within about 10% over a wide range of 
scales improving to 4.5% on the smallest scales. The lower 
panel shows a similar error bar prediction for estimates of 
the band -averaged b(k). This plot assumes r g to be fixed at 
r g — 1 over all scales. The bias parameter b(k) is even more 
tightly constrained at a few percent over a wide range. 



10 CONCLUSIONS 

In this paper we have presented the formalism for the indi- 
vidual and combined Fisher information analysis for galaxy 
redshift and velocity field surveys. This analysis allows us 
to optimise both surveys to maximise the information con- 
tent for cosmological parameters, providing an estimate of 
the uncertainty on the measurement of the matter and ve- 
locity power spectra and the set of cosmological parameters, 
(A g , r, P, u>\,, r g ). For both 6dF redshift and velocity surveys 
we find the optimal design to be as wide as possible - a result 
which was previously well known. In the case of the velocity 
survey we find the best design to be as well sampled and as 
accurate as possible and in the redshift survey we find an 
optimal sampling of around 70%. We expect to be able to 
constrain A g , V and f3 to around 2 — 3% from the redshift 




Figure 14. Constraining scale dependence: Ab(fc) and Ar 9 (fc). 



survey. From the velocity survey A v can be constrained to 
5% but a joint constraint of A v and F will have marginalised 
uncertainties of 25%. 

We find that the major benefits of 6dF are found when 
the velocity and redshift surveys are combined and when 
we wish to jointly constrain the parameters r g and j3. The 
parameters' degeneracy is broken when the power spectra 
are combined and the parameter r g can be measured much 
more accurately than in any of the above surveys with just 
redshifts. Finally, the scale dependence of r g and b can be 
measured with the combined data set - which at least will 
give credance to some of the assumptions commonly made 
about biassing. Clearly the great benefit of peculiar velocity 
information is that it tells us about the underlying mass and 
by combining this information with galaxy redshifts we can 
learn much about the relationship between luminous and 
dark matter. 
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APPENDIX A: JOINT ESTIMATES OF THE GALAXY CLUSTERING AMPLITUDE, A g , AND THE 
REDSHIFT-SPACE DISTORTION PARAMETER, /3 

A.l The redshift-space Fisher Matrix 

In this Appendix we derive some exact formulae for estimating the uncertainty on a joint measurement of the amplitude of 
galaxy clustering, A g , and the redshift distortion parameter, f3, in the limit of negligible shot-noise and complete correlation 
between galaxies and matter, r g — 1, and all other parameters are known. We model the survey volume by an effective 
constant volume, V . The Fisher matrix in this limit is 

Tii = \ V j (C^ 9i ln P " ( *' ^ dj ln M) (47) 

were i and j take the values A g and (3. 

The galaxy redshifted power spectrum with r g = 1 is given by (c.f. equation (18)) 

P° g (k) = D 2 (ka v ^(l + 2^(3 + f i 4 p 2 )b 2 P mm (k). (48) 

Differentiating and angle averaging we find the Fisher matrix is 
T = 



T Ag A g T&A g 

2/A 2 g <S>(f3)/A g f3 \ k^V 
$(P)/A g Q(f3)/P 2 ) for* 



(49) 



where 



.,„,-, = r\(^\ =i u.^i\ (51) 



i 



Inverting the Fisher matrix we find the parameter covariance matrix is 



T- 1 



((AA g ) 2 ) (A(3AA g ) 
(A(3AA g ) ((A/3) 2 ) 

A 2 g e(/3) -A g p${p) \ Gtv 



-A flj 8*(/3) 2(3 2 J fc3y*(/?) 

(0.01) 2 A 2 -(0.02) 2 A g f3 \ ( k \ " / V 



(0.02) 2 A 9 /3 (0.06) 2 /3 2 I V 0.2hMpc-V V 1Q8 [ h^Mpc] 



(52) 



where 



*03) = A 2 9 P 2 Det ^ = 2 + ~ 2 Ml^M!) (53 ) 

and the approximate numerical results in the last line of equation (52) are for (3 = 0.5. From this we can read off the 
marginalised fractional error on the amplitude of galaxy clustering, AA g /A g = 0.01, and the distortion parameter, A/3//3 — 
0.06, in the absence of shot-noise. For a single-parameter estimate of (3 the fractional uncertainty reduces to equation (36). 

The correlation coefficient between the redshift-space distortion parameter and the amplitude of clustering in the limit 
of negligible shot-noise and all other parameters known is 

-•*■*= -imr-" m - (54) 

where the final numerical result is again assuming (3 = 0.5. 
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A. 2 The uncertainty on an estimate of the amplitude of mass clustering 

Taylor et al (2001) and Tadros et al (1999) proposed that the amplitude of mass clustering can be estimated from redshift 
surveys, if galaxies are perfectly correlated with mass r g = 1, from 

Q = (A g f3)n m °- 6 - (55) 

If the errors on A g and /3 are similar, then this combination has the advantage that it can be estimated more accurately than 
either of its parts because it marginalises along the longest eigenvalue of the parameter covariance matrix. We can estimate 
the uncertainty from this combination as 



AQ / {(A4) 2 ) {JAm | g (AA g A/3) 



6tt 2 / 9(/3) - 2$(/3) + 2 
k 3 V \ *(/?) 

For f3 = 0.5 this reduces to 



(56) 



Q ' V°- 2/lM P c_ V v i08 ( ft_1M P c ) 3 

Hence although the projected uncertainties on A g and /3 are expected to be 1% and 6%, the fractional uncertainty on the 
more physical parameter Q, can be determined to 5% from redshift surveys, under the assumptions stated above. This is a 
better constraint than j3 but not A g . Because the error bars in A g and j3 are so different, the major axis of the error ellipse 
does not lie at 45 degrees to the parameter axis and so the marginalisation does not allow Q to be better constrained than 
both of the independent parameters. 



APPENDIX B: THE BIVARIATE FISHER MATRIX FOR TWO CORRELATED GAUSSIAN FIELDS 

In this Appendix we explicitly calculate the Fisher matrix for two correlated Gaussian fields. If we label the fields X and Y 
the Fisher matrix may be written in the general form: 



Ta= f d 3 fcd 3 r^a i C*x[Cov(C*xC* Y )]" 1 a,C*Y 



(58) 



(Zaldarriaga, Spergel and Seljak 1997). In our case, X and Y can denote gg, u'u' and gu' . We can rewrite equation (58) in a 
way that isolates an effective volume: 



Ta = J tfk^diCxdjCY J d 3 r[Cov(C x CY)] _1 



(59) 



The diagonal components of the spectral covariance matrix are: 

Cov(C gg (k,r)C gg (k,r)) = 2[P gg (k) + N gg (r)] 2 , (60) 
Cov(CWCW) = 2[P u , u ,(k)+N u , u ,(r)} 2 , (61) 
Cov(C u , g (k,r)C u , g (k,r)) = [(P u , g (k)) 2 + (P u , u ,(k) + N u , u ,(r))(P gg (k) + N gg (r))] (62) 

with off-diagonal parts 

Cov(C u , u ,C gg ) = 2(C gu ,f, (63) 

Cov(C„/ u iC fl „0 = 2C gu ,(C u , u ,+N u , u ,), (64) 

Cov(C gg C gu ,) = 2C gu ,(C gg + N gg ). (65) 

The full expression for the bivariate Fisher matrix is then: 
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where 



d A k 



— / ^2^)3 J d r (C u i u i diCggdjCgg + 2C gg C u i u i diC gu i djC gu i -\- C gg diC u i u i djC u i u I 

— 2C u l u >Cg U l {diCgu'djCgg + diCggdjCg U '] 



+Cg U i[diC u i u /djCgg + 2diC'g U idjC' gu i + diCggdjC u i u /])/(2(CggC u / u i — Cg U /) 2 ) (66) 



Cgg = Pgg(k) + N gg (r), 



Cg U , = Pu'g(k) =r gV / Pgg(k)P u , u ,(k). (67) 
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